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NOTATION 


The  final  forms  of  the  equations  in  the  analysis  used  only  dimension- 
less coefficients  and  ratios.  However,  during  the  derivation  the  U.S. 
Customary  System  of  Units  was  used  as  an  anchor  for  maintaining  consis- 
tency of  units.  Equivalent  values  in  the  international  System  of  Units 
(SI)  are  indicated  parenthetically. 

A, B Wave  drag  curve  fit  coefficients 

Aj.  Channel  entrance  area,  ft2  (m2) 

Aj  Fan  disk  area,  ft2  (m2) 

A Area  of  sidewall  gap,  ft2  (m2) 
s 

A^  Area  of  trailing  edge  gap,  ft2  (m2) 

A Submerged  or  emersed  area  of  sidewall,  ft2  (m2) 
w 

C „ Coefficient  of  vehicle  form  drag,  d /*5p  V 2S 
DP  f a “ 

C,  Friction  coefficient  for  water,  d /*sp  V 2A 
f w w 00  w 

C.  Friction  coefficient  for  air,  d /^p  V 2S 
fa  a a °° 

C,  Lift  coefficient,  L/4p  V 2S 

L a 00 

C Pressure  coefficient  under  the  wing,  ^ 

p *77* 

c Chord  length,  ft  (m) 

D Vehicle  drag,  lb  (N) 

D Wave  drag,  lb  (N) 

w 

d Air  drag,  lb(N) 

a 

dj  Form  drag,  lb  (N) 

d Hydrodynamic  drag,  lb  (N) 

w 

F Froude  number,  V/  /gcT 

f(h)  Cruise  height  function 


v 


i 


2 2 I 

g Gravitational  constant,  ft/sec  (m/sec  ) 

H Height  of  wing  above  water,  ft  (m) 

1 

h Cruise  height,  ft  (m) 

h Wave  height,  ft  (m) 

w 

L Lift,  lb  (N) 

j 

i.  Characteristic  length  - chord,  ft  (m)  • 

m Wave  drag  hump  factor 


n .-J 


ra,.  Forward  mass  flow  rate  - A.p  V. 
f raj 

m.  Mass  flow  rate  of  fans  = A o V, 

J j a j 

m Mass  flow  rate  from  under  wing  trailing  edge  = A p V 
te  te  a j 

ra  Mass  flow  rate  from  under  sidewalls  **  A p V 
s s a cr 

n Wave  drag  curve  slope 

2 2 

p Local  static  pressure,  lb/ft  (N/m  ) 

2 2 

pro  Free  stream  static  pressure,  lb/ft  (N/m  ) 

2 2 

S Wing  planform  area,  ft  (ra  ) 


Thrust,  lb  (N) 

2 

Critical  velocity  * V - V , ft/sec  (m/sec) 

J 3 

Jet  velocity,  ft/sec  (m/sec) 

Velocity  under  the  wing,  ft/sec  (m/sec) 


V Free  stream  velocity,  ft/sec  (m/sec) 

oo 


W Vehicle  weight,  lb  (N) 

3/2 

WT  Vehicle  density  = (W/S  ) (1/p  g) 

w 

p Density  of  air,  lb/ft^  (N/m^) 

3 3 

p Density  of  sea  water,  lb/ft  (N/m  ) 
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ABSTRACT 


Recent  efforts  in  theoretical  analysis  and  experimental  observations 
have  moved  the  concept  of  power-augmented  flight  in  ground  effect  toward 
practicality.  With  the  addition  of  lift  and  drag  due  to  external  airflow, 
end  plate  leakage,  water  skin  friction,  wave  drag  effects,  and  wave 
clearance  constraints,  a flow  model  can  be  made  capable  of  effective  com- 
parison of  different  vehicle  configurations  in  various  operating  environ- 
ments. The  analysis  shows  that  there  are  wave  drag  and  wave  clearance 
related  limits  to  many  aspects  of  the  vehicle  configuration.  The  most 
important  of  these  are  cruising  height  and  aspect  ratio,  in  that  they 
have  a very  large  effect  on  transport  efficiency.  Unfortunately,  the 
low  flying  high  aspect  ratio  cases  are  ruled  out  du&vtcuV&rtr  impact 
problems.  Testing  of  a point  design  vehicle,  arrived  at  through  the  use 
of  an  analysis  such  as  the  one  of  this  report,  needs  to  be  performed  as 
final  verification  of  the  accuracy  of  the  design  procedure. 


ADMINISTRATIVE  INFORMATION 


This  report  supports  the  Advanced  Naval  Vehicle  Concept  Evaluation 
(ANVCE)  conducted  by  the  ANVCE  Project  Office  (NOP-96V) . The  Aviation 
and  Surface  Effects  Department  (ASED)  of  the  David  W.  Taylor  Naval  Ship 
Research  and  Development  Center  (DTNSRDC)  did  the  work.  The  Naval  Air 
Development  Center  (NADC)  funded  the  effort  under  Task  Area  SSH15002 
with  Project  Order  N62269/77/PO/O0568.  The  DTNSRDC  Work  Unit  Number 
was  1612-009. 


INTRODUCTION 


Inherent  in  the  development  of  any  class  of  vehicle  is  the  need  to 
predict  their  performance  characteristics.  Only  recently  has  there  been 
the  necessary  development  of  theory  (Ref  1-2)  and  substantiating  experi- 
ments on  power-augmented  ram  wings  to  allow  this.*  With  the  addition  of 
appropriate  wave-drag  factors  derived  from  References  6 and  7 and  con- 
straints on  such  things  as  wave  clearance  and  minimum  C^  requirements, 
a vehicle  analysis  program  has  been  generated. 


*Reported  informally  by  F.  Krause  ("Static  Performance  of  a Power-Augmented 
Ram  Wing,"  DTNSRDC  ASED  TM-16-76-74 , May  1976;  and  "Parametric  Investigation 
of  a Power-Augmented  Ram  Wing  Over  Water,"  DTNSRDC  ASED  TM-16-76-95,  Oct 
1976). 
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The  computer  program  described  in  this  report  has  the  capability  of 
investigating  a large  number  of  vehicles  of  varying  configurations 
operating  in  a variety  of  environments.  It  can  also  investigate  the 
off-design  performance  of  any  configuration  that  proves  promising. 

THEORY 


FLOW  MODEL 

External  flow  effects,  end-plate  leakage,  end-plate<wetting~cfrag, 

propulsor  scrubbing  drag,  and  pressure  patch  hydrodynamic  wave  drag  are 
added  to  the  2-D  static  Power-Augmented  Ram  Wing  (PAR)  flow  models  of 
References  1 and  2.  These  earlier  models  were  corroborated  by  2-D 
experiments  for  geometries  giving  high  performance.  Also,  the  basic  mass, 
momentum,  and  energy  balances  on  which  the  static  theory  rests  are  as 
valid  in  3-D  as  for  the  2-D  case  used  in  the  original  derivation.  This 
3-D  validity  has  been  corroborated  experimentally.* 

The  external  flow  effects  amount  to  additional  lift  and  drag  propor- 
tional to  the  free-stream  dynamic  pressure.  These  forces  are  calculated 
by  assigning  constant  external  lift  and  drag  coefficients.  Presumably, 
these  coefficients  can  be  bounded  within  suitably  narrow  ranges  allowing 
concentration  on  the  PAR  features  of  the  design. 

End-plate  leakage  detracts  from  performance  by  allowing  the  air  to 
escape  from  the  cushion  in  a direction  not  parallel  to  the  free-stream 
flow.  This  is  basically  the  same  mechanism  as  thrust  loss  due  to  diver- 
gence in  conical  rocket  nozzles.  As  the  air  passes  under  the  end  plates, 
it  is  assumed  to  accelerate  laterally  only  and  to  retain  the  same  stream- 
wise  component  of  velocity  it  had  under  the  wing.  This  model  essentially 
requires  separation  and  the  concomitant  shedding  of  a vortex  sheet  at  the 
bottom  of  the  end-plate.  Since  both  of  these  phenomena  have  been  ob- 
served experimentally  on  unaugmented  ram  wings  (Ref  3),  there  is  no 
reason  to  suppose  that  they  would  not  occur  here. 

*Ibid.,  p.2. 
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The  drag  of  the  end-plates  cutting  through  the  tops  of  water  waves 
is  calculated  by  assigning  a turbulent  skin  friction  coefficient  and 
multiplying  by  the  wetted  area  and  the  dynamic  pressure  of  the  water. 

Certain  empirical  data  (Ref  4)  and  test  data  (Ref  5)  suggest  that  wetting 
drag  could  be  much  less  at  high  speeds.  End-plates  with  higher  drag  char- 
acteristics could  also  be  constructed.  These  refinements  however,  are 
ignored.  End-plate  hydrodynamic  lift  is  also  neglected. 

Wave  drag  is  calculated  from  a curve  fit  to  the  Newman  and  Poole  (Ref  6) 
analytical  results.  Subsequent  calculations  by  Doctors  (Ref  7)  indicated 
'that'  Tower" wave  d fUg" T hroUgh**4ma/p*  is -i'iiteiy-d-iie  4o-‘-t-he'  giradual^piSessUre  . 
changes  at  the  leading  and  trailing  edges  in  the  real  case,  compared  to 
the  abrupt  pressure  changes  assumed  by  Newman  and  Poole.  This  refinement, 
however,  is  inconsequential  at  the  high  cruise  speeds  of  the  PAR  wing-in- 
ground effect  (WIG)  vehi'le. 

In  summary,  the  momentum  fluxes  described  in  Figure  1 are  summed  to 
calculate  the  net  thrust  on  the  system.  The  various  momentum  fluxes  at 
the- leading  edge  are  related  by  the  previously  developed  2-D  solutions. 

The  flow  under  the  end-plates  and  trailing  edge  flap  is  modeled  as  a 
simple  one-dimensional  flow.  In  equilibrium,  the  resulting  lift  is  equal 
to  the  constant  pressure  acting  on  the  bottom  of  the  wing  times  the  wing 
area,  plus  the  aerodynamic  lift  due  to  the  external  flow. 

The  power  required  by  the  system  is  the  "air  power"  of  the  propulsor, 
defined  as  the  mechanical  energy  flux  through  the  propulsor  exit  minus 
the  mechanical  energy  flux  into  the  propulsor  inlet. 


ANALYSIS  AND  DERIVATIONS 


The  theory  used  in  this  analysis  was  based  on  the  work  of  Gallington 
and  Chaplin  (Ref  1 and  2).  The  resulting  equations  of  primary  importance 


to  this  study  are  (Figure  1): 
T = -m,V.  + 
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and  then  substituted  into  Equation  (4)  to  get  : 
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The  inverse  of  Equation  (9)  is  the  transport  efficiency,  which  is 

the  primary  point  of  comparison  between  one  vehicle  geometry  and  another. 

A unique  feature  of  the  computation  procedure  is  the  selection  of  the 

depth  of  end-plate  clipping  which  results  in  the  minimum  effective  drag. 

As  the  end -plate  runs  deeper  in  the  water,  air  leakage  under  the  en-i 

plate  is  reduced  (a  favorable  effect)  while  the  water  drag  increases  (an 

unfavorable  effect).  Refer  to  Figure  2 and  imagine  the  end  plate  being 

lifted  by  an  amount  dh  in  a wave  of  amplitude  one.  There  is  a reduction 

2 -1 

of  water  drag  proportional  to  4C^  hp  V cos  (h)dh,  and  an  increase  in  air 
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The  minimum  effec- 


drag  (or  thrust  reduction)  of  2p  V (1-V  )cos  ^(h)dh. 

a cr  u 

tive  drag  resulting  from  this  trade-off  occurs  when 


fV  ) 


f(h)  = C, 

fw  p 


cos_1  (h)  - /C~ [ 1- (1-C  )5i|cos_1(-h) 
P P 


(ID 


This  transcendental  equation  is  solved  for  h by  a Newton-Raphson  iteration 
in  the  design  procedure  of  the  program  (System  I) . In  the  analysis  pro- 
cedures (Systems  II  and  III),  h is  specified  and  the  Newton-Raphson 
procedure  is  not  used. 

The  area  available  for  side  flow  under  the  end  plates  A^  and  the 

total  wetted  area  of  the  end  plates  A are  calculated  from 

v w 

A 

/ = ^ f"”0053  1 + 2 sin  [cos  ^ 

and 

A 

sT  = f tI  sin  fcos  1 " 2 h cos  ^ 


(12) 


(13) 


The  average  size  of  the  leakage  areas  under  the  side  plate  and  the 
amount  of  water-wetted  end  plate  surface  area,  as  functions  of  the  ele- 
vation of  the  bottom  of  the  end  plate  with  respect  to  mean  sea  level, 

are  calculated  by  representing  the  design  sea  state  with  sinusoidal 
waves  corresponding  to  the  significant  waves.  The  adequacy  of  this 
assumption  may  be  assessed  by  computing,  for  the  "statistically  correct" 
sea  state,  the  probability  (averaged  over  time  or  distance)  of  encoun- 
tering water  at  each  elevation  as  a function  of  elevation,  and  comparing 
with  a similar  calculation  for  the  assumed  sinusoidal  waves.  Such  a com- 
parison is  shown  in  Figure  3 for  upper  sea  state  5.  Roth  models  show 


the  occurrence  of  air  and  water  to  be  equally  probable  at  mean  sea  level. 
The  sinusoidal  model  fits  well  at  large  and  small  elevations  and  gives 
exactly  the  correct  slope  passing  through  mean  sea  level.  Minor 
differences  occur  between  2 and  3 feet,  and  they  are  totally  negligible 
in  the  calculation  of  drag-especially  considering  the  other  approximations 
that  had  to  be  made. 


DESIGN  PROCEDURE 


SYSTEM  I 

This  system  is  intended  to  examine  the  capabilities  of  a wide  range 
of  vehicle  configurations  and  operating  parameters  in  a given  environment 
That  environment  is  definec,  ov  th"  relative  wave  height,  wave  clearance 
requirements  (i.e.,  1/3  highest,  1/1000  highest,  etc.),  and  vehicle  den- 
sity, which  is  related  to  payload-range  requirements. 


For  a given  aspect  ratio,  the  pressure  coefficient  under  the  wing 

'C  ) is  varied  from  0.1  to  0.9,  while  the  ratio  V /V.  is  varied  over  the 
P 00  .1 

same  range.  Variations  in  C are  related  to  changes  in  the  trailing  edge 

P 

gap  (and  therefore  the  flap  angle),  while  the  changes  in  V^/V^  correspond 
to  different  engine  operating  conditions.  At  each  pair  of  and  V^/V 
points,  the  optimum  performance  (or  minimum  drag)  case  is  calculated  for 
comparison  with  the  other  operating  cases  for  the  same  aspect  ratio.  Not 
all  operating  cases  are  available  for  comparison,  however,  because  certain 
constraints  (discussed  below)  might  not  be  satisfied.  Of  the  remaining 
operating  cases,  the  one  with  the  best  transport  efficiency  is  chosen 
as  the  best  case  for  that  aspect  ratio  (under  the  constraints  of  the 
specified  environment.) 


Some  computed  operating  points  are  impractical  because  of  two  possible 
limitations:  First,  it  is  difficult  to  stabilize  an  unaugmented  wing  in 

ground  effect  at  very  low  lift  coefficients  because  the  lift  becomes 
more  weakly  dependent  on  ground  clearance.  The  PAR-WIG  eventually  becomes 
unaugmented  at  the  higher  speeds  where  the  jet  velocity  and  the  free- 
stream  velocity  are  not  much  different.  Therefore,  it  seems  appropriate 
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to  specify  a lower  limit  to  the  permissible  lift  coefficient.  This  lift 
coefficient  limit  can  also  be  used  as  a correctly-scaled  structural 
integrity  speed  limit.  Second,  while  the  drag  is  calculated  using  an 
average  sinusoidal  wave,  there  are  larger  waves  in  the  design  sea  state, 
and  the  bottom  of  the  wing  must  clear  these  larger  waves.  Therefore,  we 
reject  designs  in  which  the  filled  duct  solution  indicates  a wing  height 
less  than  a certain  factor  times  the  design  wave  height.  Typically,  the 
drag  calculation  is  based  on  sinusoidal  waves  of  amplitude  equal  to  the 
significant  wave  height  of  the  design  sea  state,  and  the  wing  bottom  is 
required  to  clear  a wave  of  twice  the  significant  wave  height  (or 
the  largest  wave  expected  in  1000  wave  encounters).  Alternatively,  the 
drag  can  be  based  on  the  average  wave  height,  and  the  wing  bottom  can 
be  required  to  clear  a wave  of  three  times  the  average  wave  (again  the 
largest  wave  expected  in  1000  wave  encounters). 

This  process  is  carried  out  for  a range  of  aspect  ratios.  It  is  there- 
fore possible  to  compare  a vast  number  of  vehicles  operating  in  a variety 
of  environments.  The  configurations  that  appear  to  be  the  most  promising 
are  then  examined  in  SYSTEM  II  and  SYSTEM  III  for  their  off-design  per- 
formance . 

OFF-DESIGN  PERFORMANCE 

Having  selected  the  major  parameters,  using  the  design  procedure 
described  above,  the  next  step  is  to  predict  how  that  particular  craft 
will  perform  at  off-design  conditions.  There  are  two  off-design  condi- 
tions of  primary  interest.  First,  one  may  want  to  accelerate  as  fast 
as  possible  from  rest  to  the  cruise  condition  in  the  design  (or  lower) 
sea  state  by  keeping  the  side  plates  at  a depth  which  results  in  minimum 
effective  drag  (System  II).  The  second  important  off-design  condition 
is  operation  at  an  arbitrary  height  either  to  decelerate  by  dropping 
the  end  plates  in  the  water  or  to  clear  a given  wave  condition  with  the 
bottom  of  the  wing  (System  III).  A corollary  possibility  is  to  predict 
the  perfqrmance  of  a given  configuration  over  smooth  water  at  a given 
height  (or  end  plate  depth). 


SYSTEM  II  - EXCESS  THRUST  WITH  END  PLATES  AT  OPTIMUM  DEPTH 

For  each  given  C and  V /V.,  the  calculation  of  the  ideal  jet  area 
P * J 

for  the  filled  duct  proceeds  exactly  as  before.  For  off-design  condi- 
tions, however,  one  cannot  change  the  engine  exit  area  at  will.  The  dif- 
ference between  the  ideal  jet  area  and  the  Installed  jet  area  represents 
an  excess  thrust  or  drag,  depending  on  which  is  larger.  There  is  excess 
thrust  if  the  installed  jet  area  is  larger  than  the  optimum  jet  area. 

For  the  purpose  of  this  calculation,  we  assume  that  the  momentum  flux 
associated  with  the  excess  Installed  jet  area  compared  to  the  jet  area 
required  by  the  filled  duct  solution  appears  as  excess  thrust.  This  ex- 
cess thrust  is  then  used  to  compute  the  acceleration.  There  is  an  under- 
lying assumption  that  the  theoretical  filled  duct  performance  is  achieved 
for  the  calculated  (required)  jet  area.  In  Figure  2,  a comparison  is 
made  between  the  design  flow  model  and  the  off-design  flow  model. 

SYSTEM  III  - PERFORMANCE  AT  ARBITRARY  HEIGHT 

In  System  III,  the  acceleration  or  deceleration  of  the  craft  operating 
at  a height  different  from  that  which  would  result  in  minimum  drag  is 
predicted.  The  calculation  procedure  is  exactly  as  in  System  II,  except 
that  the  end-plate  wave  cutting  depth  is  not  optimized  for  minimum  drag; 
Instead,  the  wave  cutting  depth  is  set  at  some  specific  value.  Again, 
the  required  jet  area  for  the  filled  duct  solution  is  computed.  If 
the  installed  jet  area  exceeds  the  calculated  jet  area,  the  excess  thrust 
is  calculated  exactly  as  in  System  II  (see  Figure  U) . 

RESULTS 

To  generate  output  suitable  for  discussion,  we  arbitrarily  assign 
the  following  constant  values: 


Lift  Coefficient  of  External  Flow 

CLU 

= 0.3 

Turbulent 

Air  Skin  Friction 

Cfa 

= 0.003 

Turbulent 

Water  Skin  Friction 

Cfw 

- 0.0003 
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Drag  Coefficient  of 

External  Flow 

CDP 

- 0.04 

Minimum  Cruise  Lift 

Coefficient 

^Lmin 

- 1.5 

Minimum  Wing  Clearance 

H 

*=  2 (h 

h 

w 

f> 

H 

= 3 (h 

h 

w 

f> 

l/y  used 
for  drag) 

AVE*  USed 
for  drag) 


GENERAL  TRENDS 

All  of  the  results  presented  represent  the  best  performance  conditions 
available  for  the  parameters  specified,  including  wave  clearance  and  C 

Li 

limits . 

An  increase  in  aspect  ratio  results  in  a general  increase  in  the 
transport  efficiency  regardless  of  changes  in  other  parameters  (see 
Figures  5,  7 and  9).  Conversely,  as  aspect  ratio  increases,  the  Froude 
number  based  on  /~S  decreases  to  a greater  or  lesser  degree;  see  Figures 
6,  8,  and  10.  Again,  this  general  trend  is  evident  regardless  of  the 
variation  of  other  parameters  such  as  wave  clearance,  vehicle  density, 
and  relative  wave  height.  The  wave  clearance  constraints,  which  are 


for  the  irregular  (almost  quantum)  jumps  in  these  curves.  This  is 
because  an  increase  in  aspect  ratio  results  in  a decrease  in  the  cruise 
height,  since  the  frontal  area  (defined  by  span  x cruise  height)  is 
held  constant.  At  some  point,  this  cruise  height  will  be  less  than  the 
specified  wave  clearance,  and  it  will  be  ruled  out  as  a possible  operating 
condition. 

As  can  be  seen  in  Figure  5,  the  transport  efficiency  increases  with 
decreasing  relative  wave  height.  This  is  because  the  WIG  can  cruise 
closer  to  the  surface  with  less  interference  from  wave  impact.  Figure  6, 
however,  shows  that  the  Froude  number  decreases  with  decreasing  relative 
wave  height.  The  reduction  in  best  Froude  numbers  (based  on  rS)  in  smaller 
waves  is  probably  due  to  a reduction  In  thrust  loss  under  the  end  plates, 
allowing  a lower  efficient  cruise  speed  as  compared  to  that  in  higher 
waves. 
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Vehicle  density  has  virtually  no  effect  on  transport  efficiency 
(Figure  7),  but  it  has  a clear  effect  on  the  Froude  number.  Figure  8 
shows  an  increase  of  20  to  30-percent  in  the  Froude  number  when  the  density 
of  the  vehicle  is  doubled. 

Only  a 10-  to  25-percent  increase  in  transport  efficiency  was 
achieved  when  the  wave  clearance  constraint  was  relaxed  by  50-percent 
(Figure  9).  The  effect  of  that  same  change  in  the  clearance  constraint 
on  the  Froude  number  was  mixed  and  showed  some  dependence  on  aspect  ratio 
(Figure  10).  Above  an  aspect  ratio  of  2.0,  the  lower  wave  clearance 
allows  a lower  Fourde  number  (e.g.,  a lower  cruise  velocity  and  cruise 
height),  and  therefore  a higher  transport  efficiency.  For  aspect  ratios 
less  than  2.0,  the  reverse  appears  to  be  the  case;  however,  the  trend  is 
not  sustained  enough  to  allow  great  confidence  in  that  conclusion. 


EXAMPLE  OF  OFF-DESIGN  PERFORMANCE  CALCULATION 

Figure  11  illustrates  the  calculations  of  off-design  performance  of 
a given  power-augmented  WIG.  The  performance  of  the  PAR  feature  at  for- 
ward speed  is  represented  by  the  superposition  of  lines  of  constant  C^, 

constant  V /V,,  and  constant  Froude  number.  Tnese  three  families  of 
” J 

lines  can  be  constructed  from  the  computer  output  of  System  II  for  the 
minimum  drag  assumption  or  from  System  III  for  the  constant  height 
assumption.  A "vehicle  characteristic"  can  also  be  constructed  passing 
through  all  trim  points  (T  = D,  no  longitudinal  acceleration)  in  any 
sea  state.  The  vehicle  characteristic  for  the  design  sea  state  passes 
through  the  design  point  D.  The  off-design  vehicle  characteristics  are 
obtained  by  running  System  II  or  III  for  the  given  vehicle  in  the  new 
sea  state.  The  propulsion  system  has  a characteristic  variation  of  jet 
velocity  with  forward  speed  (i.e.,  Froude  No.);  this  propulsion 
characteristic  for  the  design  case  engine  operating  at  design  power 
passes  through  the  design  point.  High  power  settings  or  low  power 
settings  are  represented  by  other,  roughly  parallel,  characteristics. 

To  illustrate  the  computation  of  off-design  performance,  consider 


* 
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the  possibility  of  variations  away  from  the  design  point  D in  Figure  11. 

As  power  is  increased  at  constant  speed,  the  operating  point  moves  along 

a constant  Froude  number  line  to  the  new  power  characteristic,  i.e.,  point 

A.  Note  that  V.  increases  and  C decreases  in  the  process  while  V re- 
J P 

mains  constant.  The  vehicle  then  accelerates,  according  to  the  values 
given  in  the  computer  output  as  a function  of  and  V^/V^,  along  the 
new  power  characterist ic  to  point  B where  the  speed  (Froude  No.)  is  higher. 
A deceleration  resulting  from  a power  reduction  occurs  in  a converse  way 
moving  first  from  D to  C and  then  to  E. 

Another  possibility  is  operation  at  other  than  the  design  sea  state. 

If  the  existing  sea  state  is  less  than  the  design  case  and  constant  power 
is  maintained,  the  craft  accelerates  (according  to  the  values  from  the 
output  of  System  II  or  III  for  the  reduced  sea  state)  along  the  propulsion 
characteristic  to  point  F.  Conversely,  if  the  sea  state  increases,  the 
vehicle  decelerates  to  G along  the  propulsion  characteristic  by  a similar 
process.  Another  possibility  is  to  reduce  power  and  maintain  speed  in 
the  smaller  waves  thus  moving  to  point  C.  The  transport  efficiency  at 
F is  generally  better  than  at  C;  this  suggests  that  the  higher  speed, 
rather  than  the  reduced  power  cruise,  is  preferred  in  the  low  sea  state. 

COMPUTER  PROGRAM 

The  performance  analysis  program  consists  of  a main  program  segment, 
three  primary  subroutines,  and  six  subordinate  subroutines.  The  main 
program  does  little  more  than  read  in  the  data  and  control  the  use  of 
the  three  primary  subroutines  based  on  a set  of  commands  (also  input 
as  data) . 

Primary  subroutine  one  (SYSTEM  I in  Figure  12)  assumes  the  filled 
duct  solution  of  the  momentum  flux  model  and  unrestricted  cruise  height 

operation.  The  subroutine  then  calculates  the  Froude  number,  cruise 

« 

height,  transport  efficiency  (including  wave  drag  effects),  and  checks 
for  wave  clearance  and  C.  minimums  by  calling  FROUDE,  HEIGHT,  TRANS,  and 

L 

CHECK,  respectively;  see  Figure  12.  This  is  done  for  a range  of  C^  and 


V^/Vj  values;  the  results  are  then  sent  to  "MAP"  (see  Table  1 for  sample 
of  MAP).  This  procedure  Is  carried  out  for  each  of  the  aspect  ratio  cases 
specified  by  the  Input  data  to  the  main  program  segment. 

Primary  subroutine  two  (SYSTEM  II  In  Figure  12)  uses  a predetermined 
aspect  ratio  and  Aj/S  ratio  In  conjunction  with  an  essentially  fixed  end- 
plate  size.  SYSTEM  II  compares  the  performance  characteristics  of  the 
minimum  drag  case  (for  a given  C^,  V^/V^  pair)  with  those  of  the  design 
point.  The  result  is  presented  in  terms  of  vehicle  excess  power  (non- 
diraensionalized  by  WV),  and  vehicle  acceleration  in  g's.  Subroutines 
HEIGHT,  CHECK,  and  FROUDE  are  utilized,  as  well  as  MAP  and  MAPB  which 
generates  the  "Acceleration  Map."  (See  Table  2 for  sample  of  MAPB.) 
Control  is  then  returned  to  the  main  program. 

Primary  subroutine  three  (SYSTEM  III  in  Figure  12)  operates  in  the 
same  way  as  does  SYSTEM  II.  SYSTEM  III,  however,  uses  a fixed  cruise 
height  instead  of  a calculated  height  based  on  minimum  drag.  It,  there- 
fore, calls  only  FROUDE  and  CHECK  along  with  generating  the  outputs  of 
MAP  and  MAPB. 

The  first  of  the  six  subordinate  subroutines  is  HEIGHT.  This  routine 
performs  a Newton-Raphson  iteration  of  Equation  (11)  and  returns  the 
minimum  drag  cruise  height.  This  height  is  an  end-plate  emersion  factor 
that  corresponds  to  a balance  between  hydrodynamic  drag  and  loss  of 
entrapped  air  from  under  the  end-plate. 

Subroutine  FROUDE  calculates  the  Froude  number  (as  a function  of 

chord)  for  each  pair  of  C and  V /V,  points  to  be  presented  on  the  map. 

P 00  j 

The  transport  efficiency  is  calculated  by  subroutine  TRANS  and  is 
based  on  Equations  (9)  and  (10)  of  the  "filled  duct"  momentum  flux  model 
presented  earlier. 

CHECK  performs  a simple  comparison  between  the  cruise  height  and 
the  combination  of  wave  height  and  wave  clearance  factor.  It  also  checks 
to  see  if  the  vehicle's  cruise  C^  is  above  a prescribed  minimum.  Appro- 
priate flags  are  included  on  the  efficiency  map  at  each  collocation  point. 
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Subroutine  MAP  presents  the  results  of  the  design  performance  analysis 
as  numerical  values  at  their  corresponding  collocation  points.  Each  group 
of  values  is  presented  in  the  following  order:  (from  top  to  bottom)  trans- 

port efficiency,  cruise  height  factor,  cruise  height  and  limit  flags 
(if  any),  and  Froude  number  (see  Table  1). 

Subroutine  MAPB  presents  the  results  of  the  off-design  performance 
analysis  as  numerical  values  at  their  corresponding  collocation  points. 

The  values  presented  at  each  point  are:  height  of  wing  lower  surface 

above  mean  sea  level  (divided  by  /s) , vehicle  acceleration  (in  g's), 
and  P/WV  (again  in  order  from  top  to  bottom  in  Table  2). 

CONCLUSIONS 


A computer  program  was  generated  to  evaluate  the  performance  charac- 
teristics of  power-augmented  ram  wing  vehicles.  The  program  was  based 
on  an  amalgam  of  theoretical  and  experimental  investigations  regarding 
PAR  flow  phenomena  and  wave  drag  behavior.  Off-design  behavior  of 
selected  configurations  can  also  be  evaluated  in  a variety  of  operating 
environments.  With  the  aid  of  this  computer  program  some  important 
conclusions  can  be  drawn. 

o The  effects  of  increasing  aspect  ratio  are  to  increase  the  trans- 
port efficiency  and  decrease  the  Froude  number  (based  on  /S)  . This, 
however,  results  in  a decrease  in  the  cruise  height,  and  is  subject 
to  wave  impact  constraints. 

o Vehicle  density  seems  only  to  affect  Froude  number  (in  a direct 
way),  while  leaving  the  transport  efficiency  essentially  unchanged. 

o A decrease  in  relative  wave  height  results  in  an  increase  in  the 
transport  efficiency;  cruise  height  is  reduced  and  there  is  less  flow 
from  under  the  end  plates. 

o The  results  of  the  program  analysis  seem  reliable  enough  to 
warrant  the  logical  next  step  of  any  design  procedure.  That  next  step 
should  be  to  choose  a promising  design  and  compare  the  predicted  perfor- 
mance with  that  of  a test  vehicle  (probably  a fully  articulate  model). 
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Figure  7 - Effects  of  Vehicle  Density  on  Transport 
Efficiency  as  a Function  of  Aspect  Ratio 
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Figure  11  - Off-Design  Performance 
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Figure  12  - Performance  Program  Format 


TABLE  1 - SAMPLE  TRANSPORT  EFFICIENCY  MAP 
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TABLE  2 - SAMPLE  ACCELERATION  MAP 


DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 


(1|  DTNSRDC  REPORTS,  A FORMAL  SERIES  PUBLISHING  INFORMATION  OF 
PERMANENT  TECHNICAL  VALUE,  DESIGNATED  BY  A SERIAL  REPORT  NUMBER. 

(21  DEPARTMENTAL  REPORTS,  A SEMIFORMAL  SERIES.  RECORDING  INFORMA 
TION  OF  A PRELIMINARY  OR  TEMPORARY  NATURE,  OR  OF  LIMITED  INTEREST  OR 
SIGNIFICANCE.  CARRYING  A DEPARTMENTAL  ALPHANUMERIC  IDENTIFICATION. 

(3)  TECHNICAL  MEMORANDA.  AN  INFORMAL  SERIES.  USUALLY  INTERNAL 
WORKING  PAPERS  OR  DIRECT  REPORTS  TO  SPONSORS,  NUMBERED  AS  TM  SERIES 
REPORTS;  NOT  FOR  GENERAL  DISTRIBUTION. 


